Emergence of highly degenerate excited states in frustrated magnet MgCr 2 04 
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High degeneracy in ground states leads to the generation of exotic zero-energy modes, a repre- 
sentative example of which is the formation of molecular spin liquid-like fluctuations in a frustrated 
magnet. Here we present single-crystal inelastic neutron scattering results for the frustrated mag- 
net MgCr2 04, which show that a common set of finite-energy molecular spin excitation modes is 
sustained in both the liquid-like paramagnetic phase and a magnetically ordered phase with an 
extremely complex structure. Based on this finding, we propose the concept of high degeneracy in 
excited states, which promotes local resonant elementary excitations. This concept is expected to 
have ramifications on our understanding of excitations in many complex systems, including not only 
spin but also atomic liquids, complex order systems, and amorphous systems. 

PACS numbers: 75.25.Dk, 75.40.Gb, 78.70.Nx, 63.50.-x, 36.10.-k 



The concept of elementary excitations, or quasiparti- 
cles, constitutes the basis of modern condensed matter 
physics [l|- Intricate interactions among a large number 
of particles, like in liquids, can be successfully treated as 
a collection of independent quasiparticles; examples in- 
clude magnons, phonons, and rotons arising from helium 
superfluids 0]. In addition, magnetic pseudo-monopoles, 
which have been long sought after in high-energy physics, 
were recently observed as elementary excitations in a 
spin liquid-like paramagnetic phase in a highly frustrated 
magnet called a spin ice [3j]. Thus, studying elementary 
excitations in complex systems like liquids could have 
broad implications across many fields of physics. 

Highly frustrated magnets are promising sources for 
exotic spin liquid-like states. This is because in frus- 
trated magnets, not all classical-spin pairs can be ar- 
ranged antiferromagnetically on a triangular or tetrahc- 
dral lattice, which gives rise to an inherent macroscopic 
degeneracy in ground states 0, [B| . Therefore, in a low- 
temperature paramagnetic phase, frustration suppresses 
magnetic ordering and promotes spin liquid-like fluctu- 
ations (zero-energy modes), which are accompanied by 
short-range spatial correlations that lower system en- 
tropy, as required by the third law of thermodynamics. 
A representative example of this phenomenon is the for- 
mation of molecular spin liquid-like fluctuations, where 
a spin molecule refers to a spin cluster that is spatially 
confined within a geometrical region, such as an atomic 
molecule. For example, in the spinel antifcrromagnet 
MgCr 2 4 (T N = 13 K), the magnetic ions Cr 3 + {{t 2g f ', 
spin S = 3/2) form a corner-sharing tetrahedral lattice 
called a pyrochlore lattice, as shown in Fig.QJa), and the 
paramagnetic phase exhibits antiferromagnetic spin hex- 



amers @, as shown as the first mode in Fig. QJb). The 
hexamers correspond to a characteristic neutron scatter- 
ing intensity pattern that is widely spread along the Bril- 
louin zone boundary [6| , as shown in the bottom-left cor- 
ner panel in Fig. HJb), which is reproduced by classical- 
spin Monte Carlo simulations [7[. 

Notably, there is little available information on finite- 
energy elementary excitations in this molecular spin 
liquid-like paramagnetic phase; in fact, there is no ex- 
perimental report in existence. This is probably because 
zero-energy modes are one of the most important phe- 
nomena directly resulting from frustration. Indeed, many 
studies on frustrated magnets in the paramagnetic pha se 
were focused solely on the zero-energy modes [H [8l4ll| . 

By contrast, in the magnetically ordered phase in 
several spinel antiferromagnets, where frustration was 
assumed to be relieved by a lattice distortion, molec- 
ular spin and spin-orbit resonances were recently dis- 
covered to exist as non-dispersive gapped elementary 



excitation modes [a, |lfj, LL2J. An antiferromagnetic 



phase in MgC^CU exhibits a complex tetragonal spin- 
lattice order that is essentially equivalent to that in 
ZnCr 2 4 H3 |. and the hexamer mode and a hep- 
tamer mode (the first and second modes in Fig.QJb)) are 
observed at 4.5 meV and 9.0 meV, respectively [6J. In- 
terestingly, the heptamer-type zero-energy mode is also 
observed in a paramagnetic phase in another frustrated 



pyrochlore magnet Td^T^Cv [111], and a ground state 
composed of the same structural units is observed in frus- 
trated _spinel magnet AIV2O4 with charge degree of free- 
dom [15j. Thus, the magnetically ordered phase does 
incorporate some sort of frustration effect in its excita- 
tions. 



2 



[010]* 



<111> 




(a) [100] 



1st mode 2nd mode 






-6Ji 7- ^ -4^ 
3rd mode 



(b) 



FIG. 1: (Color online) (a) Pyrochlore lattice. The bold lines 
mark examples of positions of the hexamer and heptamers 
shown in (b) among equivalent positions, (b) Snapshots of 
spin molecular models. The arrows indicate spins, which dy- 
namically fluctuate in arbitrary directions in keeping with the 
relative correlations. The total exchange interaction energies 
in the molecules are also shown, where J\ is the magnitude 
of the nearest-neighbor exchange interaction energy. 



In this study, we explored finite-energy excitations in 
the paramagnetic phase and clarified their relation with 
those in the magnetically ordered phase through a com- 
prehensive study of spin excitations distributed over wide 
momentum (Q) and energy (E) ranges both in a para- 
magnetic and magnetically ordered phase in MgCr204 
by inelastic neutron scattering. The overall study of 
spin excitations was possible because of the combination 
of a large single-crystal assembly, an advanced time-of- 
fiight spectrometer with large-solid-angle detectors, and 
sophisticated software to handle huge data sets in four- 
dimensional space. 

Inelastic neutron scattering experiments were per- 
formed using the direct geometry chopper spectrometer 



MERLIN at the ISIS (UK) spallation neutron source [lfj. 
The detector coverage is as large as —45° to 135° in 
the horizontal plane and ±30° in the vertical direction. 
The incident energy (Ei) was fixed at 50 meV with a 
chopper speed of 400 Hz. The energy resolution under 
elastic conditions was approximately 5% of Ei. Single- 
crystal rods of MgCr2C>4 were grown by a floating zone 
method. Details of the crystal growth are summarized 



in Ref. [17] . The rod was about 4 mm in diameter and 
40 mm long. Six co-aligned single crystals were fixed by 
thin aluminum plates and inserted in a closed-cycle He 
refrigerator with 4 He exchange gas. Since the spin sys- 
tem is three-dimensional, as shown in Fig. HJa), the data 
were recorded while rotating the crystal in 1° steps about 
the vertical axis, which is unlike experiments on two- 
dimensional and one-dimensional systems with a fixed 
crystal angle (w) The huge combined data sets 

were handled by the HORACE software of ISIS HJ. 

Figure [2ja) shows the E spectra of Q-integrated in- 
tensity, measured in a paramagnetic phase (22 K) and 
magnetically ordered phase (8 K). In neutron scattering, 
the magnetic scattering intensity decreases by a magnetic 
form factor while the phonon intensity increases in pro- 
portion to Q 2 with increasing Q 2l| . Since below ~ 20 
meV, the low-Q intensity is higher than the high-Q inten- 
sity, as shown in Fig. [^a) , the signals in this low-energy 
region are magnetic in origin. As shown in the left panel, 
a single quasielastic mode seems to exist above Tn. In 
contrast, as shown in the right panel, peaks around 4.5 
meV and 9.0 meV and tiny shoulders around 13.5 meV 
and 18.0 meV are observed below Tn, which are of equal 
energy interval. 

Despite the appearance of a single mode, we investi- 
gated the scattering intensity distributions in a Q space 
sliced at the four characteristic E's not only below but 
also above Tn, as shown in Figs. [2jb) and[2jc), respec- 
tively. Above Tn (Fig. EJb)), different patterns spread- 
ing over several rlu appear at 9.0 meV, 13.5 meV, and 
18.0 meV in addition to the zero-energy hexamer patterns 
at 4.5 meV, indicating the existence of other extremely 
short-range correlation modes. Further, as shown in 
Fig. Hie), all the modes are sustained below T N : the 
4.5-meV patterns are almost identical. Each 9.0-meV 
pattern above T N is regarded as a superposition of the 
4.5-meV pattern and the heptamer 9.0-meV pattern be- 
low Tn and, thus, the heptamer mode certainly exists 
both below and above Tn, and the remaining 13.5-meV 
and 18.0-meV patterns above Tn are smeared compared 
with those below Tn but are essentially the same. Here, 
the 4.5-meV and 9.0-meV patterns below Tn are consis- 
tent with the previous report @. Thus, in the paramag- 
netic phase with the zero-energy mode, we found multiple 
finite-energy short-range correlation modes, all of which 
are also in common with the magnetically ordered phase. 
The appearance of a single mode is due to superposition 
of all the modes that are broader in E compared to the 
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FIG. 2: (Color online) (a) E dependence of Q-integrated intensity measured in a paramagnetic phase (left panel) and mag- 
netically ordered phase (right panel). The integration ranges are 0</t<2, < < 2, and < I < 2 reciprocal lattice units 
(rlu) in the low-Q data and < ft < 2, 5 < fc < 7, and < I < 2 rlu in the high-Q data. The arrows indicate the magnetic 
peaks and shoulders. The errors are smaller than the symbols used. The lines are guides to the eye. (b)(c) Measured inelastic 
neutron scattering intensity distributions in (h, k, I, E) space in a paramagnetic phase (22 K) and magnetically ordered phase 
(8 K), respectively. Each set of three-abreast panels corresponds to the data for the given E measured in the hkO, hkl, and 
hk2 planes. The intensities were integrated over ranges of ±0.2 rlu and ±1.0 meV around the given I and E, respectively. The 
vertical tone indicates the linear scale intensity in arbitrary units in (b)-(d). The three first quadrants, boxed by white lines 
in the hkO planes, show areas previously reported [6j]. (d) One-to-one correspondence between calculated patterns as identified 
by the molecular models shown in Fig. [Tib). 
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ordered phase. 

Next, to extract information on the spatial correlations 
of these modes, we tried to reproduce the experimental 
Q patterns as two-body correlation functions of classical 
spins, like in Refs. 0, H[. The Watson-Freeman mag- 
netic form factor of Cr 3+ 22[ and the orientation aver- 
age over the equivalent directions were also taken into 
account. Through trial and error, heptamer variation 
models were also found for the third and fourth modes, as 
shown in Fig.[ljb), and the corresponding calculated pat- 
terns (Fig.^d)) are in agreement with those of Figs.^b) 
and[2{c). Further, the total exchange interaction energy 
evaluated from the numbers of antiferromagnetic and fer- 
romagnetic nearest-neighbor bonds is of equal interval as 
follows, which is also in agreement with experiments: the 
energy for each molecule can be evaluated to be — 6J1, 
— 4Ji, — 2 Ji, and J\ with an equal interval of 2 J%, where 
Ji denotes the magnitude of an antiferromagnetic first- 
neighbor exchange interaction energy; this first-neighbor 
exchange interaction has been reported to be predomi- 
nant over other exchange interactions in band calcula- 
tions (23j . However, it is uncertain whether these are a 
unique solution. 

To summarize the aforementioned results, the present 
experiments found a common hierarchical set of spin exci- 
tations ranging from the zero-energy mode to the finite- 
energy modes above and below Tn, which can be de- 
scribed by the hexamer and heptamer variation models. 
Now we turn to the significance of this finding. First, ex- 
istence both above and below Tn means that the molec- 
ular spin excitations are neither caused by magnetic or- 
dering nor likely to be normal spin waves arising from 
the magnetic order. In addition, the spin-lattice order 
is tetragonal [lH whereas the hexamers and hep- 
tamers are trigonal, as shown in Fig. [TJa); they are also 
different in symmetry This difference in symmetry has 
been very recently detected by ultrasound measurements 
in the paramagnetic phase: there are tetragonal spin ex- 
citations/fluctuations below ~ 3.4 meV originating from 
spin Jahn-Teller coupling with the tetragonal spin-lattice 
order and two types of trigonal spin excitations around 
~ 3.4 meV and in the higher-energy region, correspond- 
ing to the hexamer and the heptamer variations, respec- 
tively [23 | . 

Second, we focus on the finding that the finite-energy 
modes are spatially confined, as is the zero-energy mode. 
The spatially confined zero-energy mode is recognized as 
being the direct result of high degeneracy in the ground 
states, which restricts the formation of a normal wave- 
like magnetic order that is inevitably longer than the 
wavelength of the propagation vector. In analogy with 
this, the spatially confined finite-energy modes suggest 
high degeneracy in excited states, or the expanded con- 
cept of frustration in excitations. Indeed, this concept 
seems to be realized in the magnetically ordered phase, 
which exhibits a complex magnetic structure with mul- 



tiple pro pag ation vectors and 32 magnetic ions in its 
unit cell 13, [14 1 , implying that 32 spin wave modes are 
squeezed in an energy region. This situation is most 
probably a clearer example of the term "frustration ef- 
fect in excitations" mentioned in the previous report [f| . 
Similarly, in the paramagnetic phase without a magnetic 
order, since the concept of a unit cell is broken down or a 
unit cell with the larger numbers of magnetic ions can be 
defined, further higher degeneracy is expected in the spin 
excited states. Thus, one can conclude that although the 
molecular modes might have originally been normal spin 
waves, the expanded frustration transforms them beyond 
recognition. 

Last, we discuss the ramifications of this expanded 
concept. The concept is realized when a unit cell con- 
sists of many components. However, in contrast, from a 
mean-field theory for pyrochlore systems with only four 
magnetic sublattices, it was reported that only ground 
states are highly degenerated and that this is not the 
case for excited states [25j. This, in turn, suggests that 
the expanded concept can be also applied to atomic liq- 
uids and amorphous systems, which are similar to spin 
frustrated systems as follows: a spin-frustrated system 
exhibits, for example, a spin glass state because of the 
presence of trace impurities [26| . which is comparable 
to a supercooled liquid that undergoes a glass transi- 
tion upon experiencing a small impact. Further, super- 
cooled liquids and glasses exhibit similar resonance-like 
short-range vibrational modes called boson peaks. In- 
terestingly, a recent X-ray spectroscopic study revealed 
that boson peaks are originally identical to, but are then 
severely transformed from, lattice waves propagating in 
a crystal (phonons), and also suggested that the Q cor- 
relations are distributed over several rlu near the pseudo 
Brillouin zone boundary [27j . 

In summary, we studied spin excitations over a wide 
(Q, E) space in the frustrated spinel magnet MgCr 2 04 by 
single-crystal time-of-flight inelastic neutron scattering. 
A set of molecular spin excitation modes was commonly 
observed above and below Tn- This observation leads 
us to the concept of frustration in excitations, which is 
applicable both above and below Tn, and probably also in 
atomic liquids, complex ordered systems, and amorphous 
systems. 
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